Bs mesons from the lattice: Excited states 
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The energies of different angular momentum states of a heavy-light meson were measured on a lattice in 
We have now continued this study using several different lattices, quenched and unquenched, that have different 
physical lattice sizes, clover coefficients, hopping parameters and quark-gluon couplings. The heavy quark is taken 
to be infinitely heavy, whereas the light quark mass is approximately that of the strange quark. By interpolating 
in the heavy and light quark masses we can thus compare the lattice results with the Bs meson. Most interesting 
is the lowest P-wave Bs state, since it is possible that it lies below the BK threshold and hence is very narrow. 
Unfortunately, there are no experimental results on P-wave Bs mesons available at present. 

In addition to the energy spectrum, we measured earlier also vector (charge) and scalar (matter) radial distri- 
butions of the light quark in the S-wave states of a heavy- light meson on a lattice 5,. Now we are extending the 
study of radial distributions to P-wave states. 



1. Energies 

For simplicity we consider a system that con- 
sists of one infinitely heavy quark (or anti-quark) 
and one light anti-quark (quark) . The basic quan- 
tity for evaluating the energies of this heavy-light 
system on a lattice is the 2-point correlation func- 
tion 6*2 — see Fig. n]a). It is defined as 

• Pt+Trtc/Q(x,t,i + T)), (1) 

where is the heavy quark propagator and Gq 
the light anti-quark propagator. Pt is a linear 
combination of products of gauge links at time t 
along paths P and F defines the spin structure 
of the operator. The (...) means the average over 
the whole lattice. The energies are then extracted 
by fitting the C'2 with a sum of exponentials, 

C2(T) « £ c,e— (2) 

i=l 
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where the number of exponentials, A'^ax, is 2, 3, 
or 4 and the time range T < 10. C2(T) is in 
practice a 3 x 3 matrix for different fuzzings (the 
indices are omitted in Eq. |21 for clarity). This 
fuzzing enables us to extract not only a more ac- 
curate value of TOi , the energy of the ground state, 
but also an estimate of the first radially excited 
state, m2. After taking the continuum limit we 
can interpolate in the light quark mass and use 
charmed meson experimental results to interpo- 
late in the heavy quark mass to get a Bg meson. 
The energy spectrum was measured using differ- 
eirt lattices: Q3 is a 16^ x 24 quenched lattice (for 
more details see PP), DF2 is a 16^ x 24 lattice and 
DF3, DF4 are 16^ x 32 unquenched lattices. DF 
refers to dynamical fermions (more details can be 
found in [3]). Parameters of the lattices are given 
in Table n 

From some theoretical considerations ,4] it is 
expected that, for higher angular momentum 
states, the multiplets should be inverted com- 
pared with the Coulomb spectrum. If the po- 
tential is purely of the form l/r the state L-|- lies 
always higher than L— . Here L-|-(— ) means that 
the light quark spin couples to orbital angular mo- 
meirtum L giving the total j = L±l/2. (Since the 
heavy quark is infinitely heavy its spin does not 
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Table 1 

Lattice parameters. Here Q means quenched and 
DF means unqucnched (dynamical fermions). 







Csw 


K 


Q3 


5.7 


1.57 


0.14077 


DF2 


5.2 


1.76 


0.1395 


DF3 


5.2 


2.0171 


0.1350 


DF4 


5.2 


2.0171 


0.1355 




a [fm] 


mq/ms 


rom^ 


Q3 


0.179(9) 


0.83 


1.555(6) 


DF2 


0.152(8) 


1.28 


1.94(3) 


DF3 


0.110(6) 


1.12 


1.93(3) 


DF4 


0.104(5) 


0.57 


1.48(3) 



Figure 1. a) Two-point correlation function C2; 
b) Three-point correlation function C3. 





Figure 2. The energy spectrum for different angu- 
lar momentum states. The notation is explained 
in the text. 2S is the first radially excited L = 
state. Here tq = 0.525(25) fm. 
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Figure 3. Interpolation from infinite mass to 
TfiQ = rrih- The results at mq = rric are from 
charmed meson experiments. Note that the low- 
est Bs P-wave state (0+) lies below the BK 
threshold and should be very narrow. 
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play a role here.) However, in the case of QCD 
the confining potential is important for large dis- 
tances and thus L— should eventually lie higher 
than L-|- for higher angular momentum states. 
Experimentally this inversion is not seen for P- 
waves, and now the lattice measurements show 
that there is no inversion in the D-wave states ei- 
ther. In fact, the D-|- and D— states seem to be 
nearly degenerate, i.e. the spin-orbit splitting is 
very small (see Fig. Note that the spectrum 
in Fig. 121 shows an approximately linear rise in 
excitation energy with L. 

2. Radial distributions 

For evaluating the radial distributions of the 
light quark a three-point correlation function is 
needed — see Fig.^b). It is defined as 



(3) 



2S 



P+ 



D- 



D+ F+- 



We have now two light anti-quark propagators, 
Gqi and Gg2, and a probe 6(i?) at distance R 
from the static quark. We have used two probes: 
74 for the vector (charge) and 1 for the scalar 
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Figure 4. Ground state vector and scalar radial 
distributions and discretized exponential fits. 



Figure 5. Radial distributions containing the first 
radially excited state (Ref. (2) for the S-wave). 
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(matter) distribution. The radial distributions, 
x*^ (i?)'s, are then extracted by fitting the C3 with 



C3(i?,r) 



J2 Qe-"'*ix^^(i?)e-™^*^cj (4) 



— see Figs.01|31 The rn^'s and q's are from the 
best fit to C2 (Eq.|2l. Here x^^ is the ground state 
distribution and x^^, for example, is the overlap 
between the ground state and the first excited 
state. 

Fig. 0] compares the ground state vector 
(charge) distribution with the scalar (matter) dis- 



tribution. The matter distribution seems to drop 
off faster than the charge distribution. However, 
at i? = the charge and matter distributions are 
roughly equal. Also note that the P-wave matter 
distribution changes its sign. (Hence the logarith- 
mic scale in the S-wave plot but a linear scale in 
the P-wave plot in Fig. 01) In Fig. El the charge 
and matter first excited state distributions are 
compared. As expected, a node is clearly seen 
in these excited state distributions. The radial 
distributions shown here are for DF2, but the 
analogous calculations for DF3 and DF4 are in 
progress. These distributions can be now treated 
as "experimental data" that needs understanding. 
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Figure 6. Dirac equation numerical solutions, F 
and G, for a potential V = —e/R + b ■ R. Here 
e = 0.6 • fic, b = (500 MeV)^ and m = 100 MeV. 
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Figure 7. Dirac equation numerical solutions 
for the same potential as in Figure compared 
with lattice results. R is given in lattice units 
(a « 0.15 fm). Here the sign of the S-wave mat- 
ter distribution has been changed for clarity. 
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The lattice measurements, i.e. the energy spec- 
trum and the radial distributions, can be used to 
test potential models. There are several advan- 
tages for model makers: First of all, since the 
heavy quark is infinitely heavy, we have essen- 
tially a one-body problem. Secondly, on the lat- 
tice we know what we put in — one heavy quark 
and one light anti-quark — which makes the lat- 
tice system much more simple than the real world. 
We can, for example, try to use the Dirac equa- 
tion to interpret the S- and P-wave distributions. 
If the solutions of the Dirac equation are called 
G (the "large" component) and F (the "small" 
component) we get 



for the charge distribution and 
x^\R) = G{Rf - F{Rf 



(6) 



x"(i?) = G{Rf + F{Rf 



(5) 



for the matter distribution. These solutions are 
plotted in Fig. El for a potential V = —e/R + 
b ■ R, where e = 0.6 • He, b = (500 MeV)^ and 
the mass that appears in the Dirac equation is 
m = 100 MeV. The values of these parameters 
are not tuned yet, because at this stage we only 
want to show that a reasonable set of parameters 
can give a qualitative description of the lattice 
data. The qualitative agreement is surprisingly 
good — see Fig. [7\ The Dirac equation approach 



gives a natural explanation, for example, to the 
fact that the matter distribution drops off faster 
than the charge distribution. Also the change of 
sign in the P-wave matter distribution comes out 
automatically. 

3. Main conclusions 

The key points are: 

• There should be several narrow Bg states, 
in particular the O"*" that lies below the BK 
threshold. 

• The spin-orbit splitting is very small, con- 
trary to the expectation predicted by a lin- 
early rising potential that is purely scalar. 

• The radial distributions of S and P— states 
can be qualitatively understood by using a 
Dirac equation model. 

The natural next step would be to calculate the 
radial distributions for P-|-, D— , D-|- and F-wave 
states. 
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